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Abstract 

We discuss the CP-violating dimuon charge asymmetry of events pp -^ 
BBX —f ^ n X in the Standard Model. Our conclusion is that the 
asymmetry is larger than previously expected and may reach a few 
percent for the \B^, B'^j system. The analysis is also extended to the 
Two Higgs Doublet Model (Model II). 



1 Introduction 

We discuss the inclusive CP-violating dimuon charge asymmetry A = ( A^++ — 
A^__)/(A^++ + N__) in the Standard Model of events pp -^ BBX -^ p^p^X 
where the BB pair is either B^Bd or BgBs- A^++ is the number of events 
pp -^ BBX -^ p^p^X. In particular we consider "indirect CP violation (or 
CP violation in the mixing)" |]1| due to complex effects in B ^^ B mixing and 
decay. 



2 Indirect CP violation 

Let us review the standard formalism of indirect CP violation]^ but without 
making the usual approximations valid when CP violation is "small". We 
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where the matrices M and F are hermitian. The hamiltonian H itself is 
not hermitian since the B mesons do decay. The matrix elements of the 
hamiltonian are obtained from second order perturbation theory: 

Ma/3 = m6af3 + « Hsw /?) - P > ^ , (2) 

Tap = 27,Y,{oi\Hw\0{^\Hw\ (3)5{E^-m), (3) 

where Hw is the Standard Model weak interaction, Hsw is a superweak 
interaction (which we consider no further), and P denotes "principal part". 

The diagonal elements of H are assumed equal due to CPT invariance. 
The argument goes as follows: Mn is the amplitude for a B^ to remain a 
B^ . Applying CPT we obtain the amplitude for a B^ to remain a -B°, i.e. 
M22. So, if CPT invariance holds, we obtain Mn = M22 = m. Likewise, 
Fii is the probability per unit time for the decay 5° — >^ ^ ^. Applying CPT 
we obtain the probability per unit time for Yl^ ^ ^^ ■ From Equation (^) 
for this process, and changing the order of the two brackets, we obtain the 
probability per unit time for B^ -^ X]^) ^-^^ ^22- So, if CPT invariance 
holds, we obtain Fn = F22 = F. 

The solution to the equation idtp/dt = Htp with ip'^ = (i?°(t),i?°(t)) is 

B'{t) = \ {s^{t) + s4t)} i?°(0) +^YV'e'\ ^"^^^^ ~ '-^^^^ ^°^°^' 
S°(t) = i±| ■ i {s+(t) - ._(t)} 5^(0) + \ {.s^{t) + s_(t)} 5°(0), (4) 



where 



s__(t) = exp(-imt)exp(-Ft/2)exp(iAMt/2)exp(AFt/4), 

s+(t) = exp(-imt)exp(-Ft/2)exp(-2AMt/2)exp(-AFt/4), (5) 



AM - fAF _ 2 (M12 - iFi2) 



l + e 2 (M*2 - |r]:2) AM - I AF 



(6) 



The phase of (1 —£)/{! + e) is arbitrary: it can be changed by redefining the 
phase of -B°(0). Observables depend on the absolute value of (1 —£)/(! + £:), 
or equivalently on 

-Re(e) ,„, 

« = 7-772- (7) 

1 + \e\ 



For the same reason, we can multiply Mu and ri2 by a common phase-factor. 
Only the relative phase is observable: 



We introduce the notation 



AM Ar 

x = ; v = — • (9) 

r 2r ^ ^ 



Then the probability that a B decays as a B is 

f°° li^Pi k _ o |2^f 
_ Jo ll+sl 4 1-^+ -^-1 "^ 



Jo ll+il il-^+^-^-l "'^ + Jo 4l'^+ + '^-l "'' 

1 + a;2 + 2a (1 - y2) ■ 
Similarly, the probability that a B^ decays as a 5" is 



Finally, let us relate Mu and ri2 with AM and AF: 

a Im {F12/M12} 



4a2 4 + IF12/M: 
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11^ = ^^, 1 + (AF/(2AM) ) 
Mi2y l + (aAF/AM) 
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(10) 



^ l + x2-2a(l-|/2) ^ > 



:i2i 



^^ ( TT- ) = -4c^ /.\ !l,^ ^'2 > (13) 



, , -a 4AM2 + AF2 

'""(^) = 131^ ■ AMAF ' ('') 



IM12I' = ^AM^ + - "" AF^ (15) 

' ' 41-4a2 41-4a2 ' ^ > 

|ri2|' = ^^AM2 + 1 ^AF^. (16) 

' ^^' l-4a2 41-4^2 ^ > 

These equations are exact for the (5°, 5°) and (-B^, 5°) systems separately. 
As a cross check, for \a\ <^ 1 we recover the well known results AM ^ 2|Mi2| 
and A ^ -4a ^ Jm (F12/M12). 
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Figure 1: A tree level Feynman diagram contributing to the 5° meson decay 
rate Fn = F oc {B^ \ Hw \ Cl {i\ Hw I B^s)- The intermediate quarks can 
also be d, u or u. There are also lepton channels. 



3 Ti2 and M12 in the Standard Model 

From the box diagrams [0] of B^ ^^ B^ mixing and Equation (^ for the 
Standard Model, we obtain [|| 



""'" r r STT^r 



where 



x^ — llx^ + ixt 3x^ In (xi) 
^{l-xtf 2 (1 - xt) 



c^i^t)^ \::' \r - :T Z - (i^) 



Let us now consider the absorptive matrix elements given by Equation 
d^). For a = j3 Equation (^ is "Fermi's golden rule" . The tree level Feynman 
diagrams are shown in Figures |l| and ^ 

We arrive at three different predictions for the CP-violating asymmetries 
Ad and Ag in the Standard Model: 

1. If we require that intermediate quark spins and momenta match in 
the diagrams of Figure H, then we obtain negligible asymmetries as discussed 
by J. Hagelinl^. This result would be correct if quarks were observable. 

2. However quarks are not observable. Hadrons are observable. Due 
to hadronization, i.e. due to gluons, the momenta need not match at 
the quark level and we obtain the approximate results shown in Table |l|. 
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Figure 2: A tree level Feynman diagram contributing to ri2 oc 
(5° I H]Y \ {^ \ Hw I B^s)- The intermediate quarks can also be u or u. 



(All input numerical data for these calculations not specified in the Tables 
were obtained from |l[].) A similar point was emphasized by T. Altomari, L. 
Wolfenstein, and J.D. BjorkenQ] so that their "conclusion is that a reasonable 
estimate of the asymmetry lies between 10~^ and 10~^ but that neither the 
sign nor the magnitude can be reliably calculated". That conclusion is in 
agreement with Table |I|. 

3. If, due to gluon exchange in the hadronization process, it were a good 
approximation to neglect both the requirements of spin and momen- 
tum match, then we obtain approximately the results shown in Table ^j. As 
an example, suppose that the c-quark in the left hand diagram of Figure | 
has "spin up" , and the c-quark in the right hand diagram has "spin down" . 
This miss-match of spins of the quarks need not imply a miss-match at the 
hadron level, as can be seen, for example, in the case in which both quarks 
hadronize into a scalar meson. 

Which of the three predictions, if any, is correct? 

4 Discussion 

Consider diagrams of the form shown in Figure |I|. Requiring match of spins 
of the quarks (and weighting hadronic modes by a factor 3 for color) we 
obtain an inclusive branching fraction B{h -^ /iX) = 0.16 to be compared 
with the experimental value B{h — > /xX) = 0.103 ± 0.005[|l|. This is the well 
known "baffling semi-leptonic branching fraction of B mesons" H . But quark 
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Figure 3: Another Feynman diagram contributing to ri2 oc 
(5° I H-v^ I ('C I Hw I -S°)- The intermediate quarks can also be u or 
u. 



spins are not observable. Only hadrons are observable. So the hadronization 
process, i.e. gluons, enhance the decay rate Fn by a factor ;^ 2 to account 
for the observed drop in semi-leptonic branching fraction. 

Let us now consider the momentum miss-match of the spectator quarks 
in the diagrams of Figure ^ To obtain matching of momenta we need at least 
one gluon as shown in Figure ^. The amplitude of the diagram of Figure § 
is "enhanced" with respect to the amplitude of the diagram of Figure ^ by a 
factor of order 

ri2(3) _ A^(3) _ dsmBJl p p /.QN 

The numbers in parenthesis refer to the Figures, and x = s or d for the S^ or 
5° systems respectively. The factor msfB '^^is introduced (somewhat arbi- 
trarily) on dimensional grounds. The gluon propagator is Pg ~ l/(2m^|p2|)- 
We take |p2| to be of order niB/^ (see definition of p2 and x in Figure ^. 
If the outgoing quarks were free, the quark propagator P^ would be on-shell 
and would diverge. Since the outgoing quarks are confined to dimensions of 
order 1 Fermi, we replace the quark propagator by P^ ~ l/?7i^. Then the 
order of magnitude estimate of the "enhancement factor" (|19|) is 



ri2(2) ^ A.,{2) ^ 16m^m, ^ ^ 
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Table 1: Standard Model calculation of the dimuon charge asymmetries A^ 
and As, and other parameters, as a function of the angles 7 and (3 of the 
unitarity triangledl, for Kfe = 0.0036, i = 1.09, and fe = 0.18GeV. The 
following parameters define the Cabbibo-Kobayashi-Maskawa (CKM) matrix 
in the Wolfenstein parametrization|l|]: A = 0.2235, A = 0.8357, [p^ + Tff'^ = 
0.3854 and 7 = arctan(?7/p) listed above. Momenta of the s-quarks in the 
diagram of Figure are not required to match. Spins are required to 



match, rub = 4.25GeV/c2, m^ = 1.25GeV/c2, 
0.006GeV/c2, and m, = 0.115GeV/c2. 



m„ 



0.00325GeV/c^ m^ 



which is ~ 0.08 for the 5° system, and ~ 1.6 for the i?° system. 

So, which of the three predictions, if any, is correct? Prediction 1 by 
Hagelin includes a restriction on \p2\ due to a miss- match of momenta of 
the spectator quark, and also requires matching spins. We have seen that 
this momentum miss-match can be "fixed" by gluons (at no cost for the 5° 
system) . Therefore we arrive at Prediction 2 given in Table |I| which does not 
require momenta to match but still requires matching spins. We have seen 
that hadronization relaxes the need to match spins (as required by the semi- 
leptonic branching fraction of B mesons). Then we arrive at Prediction 3 
given in Table |. Since xa = 0.723±0.032|| we obtain from Table |, Ad ^ 1% 
to 3% for the allowed range of Ktfell|]. But for the 5° system the gluon in the 
diagram of Figure ^ can even enhance the asymmetry above the prediction 
given in Table ^. Note that relaxing the requirement of matching spins we 
obtain a large AF^/r as predicted in [^, leading to interesting experimental 
consequences R]. 
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Table 2: Same as Table ^ except that the spins (in addition to the momenta) 
of the intermediate quarks in Figure |^ are not required to match. For the 
first half of the Table [p^ + rj'^f'^ = 0.3854 and Kb = 0.0036. For the second 
half of the Table [p^ + n^]^/^ = 0.54 and Kfe = 0.0050. At 7 = 0° and 180°, 
Ad = and As = 0. 

5 The Two Higgs Doublet Model 

We have repeated the calculations using the Two Higgs Doublet Model 
(Model n)[§l. The results are shown in Table ^ (which requires matching 
intermediate quark spins but not momenta). Comparing with Table |l| we 
note that the predictions for the asymmetry Ad are not significantly changed 
by this extension of the Standard Model. 

6 Conclusions 



Our conclusion is that the inclusive CP-violating like-sign dimuon charge 
asymmetry in the Standard Model is larger than previously expected[Q. For 
the (-B^, -B°) system the asymmetry is probably between ^ 1% and 3%, and 
may even be enhanced by gluons above this value. Such a large asymmetry 
would be within the reach of the next run of the D0 and CDF experiments 
at the FERMILAB Tevatron collider. 
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Table 3: Same as Table |I| but for the Two Higgs Doublet Model (Model II) 
with m^ = 45GeV/c^ and tan (5 = 30, instead of the Standard Model. 
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